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1.  INTRODUCTION 


The  burning  of  a solid  propellant  is  known  to  be  a function  of  the 
pressure  at  which  the  combustion  takes  place,  and  also  a function  of 
the  time  derivative  of  the  pressure.  One  method  to  determine  that  func- 
tion is  to  burn  a sample  of  the  propellant  in  a closed  container  ("bomb") 
and  register  the  pressure-tine  curve  generated  in  the  container  by  the 
combustion.  Under  simple  assumptions  about  the  burning  process  and 
combustion  gases  the  pressure-time  data  can  be  manipulated  to  yield  the 
desired  burning  law.  Formulas  and  algorithms  which  accomplish  this  task 
have  been  published,  for  instance,  in  References  1 through  4. 

Typically  these  algorithms  are  obtained  under  the  assumption  of 
constant  temperature  throughout  the  experiment  and  simple  equations  of 
state.  (The  algorithms  may  also  contain  empirical  quantities,  as  in 
Equation  (3)  of  Reference  2.)  Because  of  their  complexity,  the  published 
algorithms  are  difficult  to  analyze  with  respect  to  the  confidence  limits 
of  the  results.  As  a consequence  all  burning  rate  data  are  published 
without  corresponding  error  estimates. 

In  the  present  pa'^er  we  shall  derive  new  formulas  for  the  treatment 
of  the  pressure-tine  data,  based  on  somewhat  relaxed  assumptions.  Under 
conditions  which  are  often  satisfied  and  which  can  be  controlled  by  the 
experiment,  these  formulas  transform  the  problem  into  a simple  curve- 
fitting task,  to  which  standard  error  analysis  can  be  applied. 

In  Section  2 of  this  paper  we  shall  list  and  discuss  the  assumptions 
on  which  the  analysis  of  the  data  is  based.  A general  outline  of  the 
method  is  given  infection  S. 

The  ideas  presented  in  Section  3 are  applied  in  Section  4 to  gases 
with  Noble-Abel  equations  of  state  and  constant  specific  heats.  That 
equation  of  state  is  frequently  used  for  interior  ballistics  calculations 
and  often  the  only  equation  available  for  a given  propellant.  The 
theory  of  Section  3 yields  for  Noble-Abel  gases  fairly  simple  formulas. 

A further  simplification  is  obtained  if  the  ratio  of  specific  heats 
of  the  combustion  gases  is  equal  to  that  of  the  inert  gas  which  fills 

^Wallace,  N.,  "New  Formulas  for  Rapid  Calculation  of  Linear  Burning  Rates 
of  Solid  Propellants",  Picatinny  Arsenal  Technical  Report  2488,  April  19S8. 

*’Fayon,  A.M.  and  Goldstein,  F.B.,  "Evaluation  of  Solid  Propellant  Ballistic 
Properties  by  Constant  Volume  Burning",  Combustion  and  Flame  23  (1966). 

t 

Osborn,  J.R.,  "Evaluation  of  Solid  Propellant  Ballistic  Properties", 
Combustion  and  Flame  20,  193  (1973) . 

^Krier,  H.,  Shimpi,  S.A.,  Adams,  M.J.,  "Interior  Ballistics  Predictions 
Using  Data  from  Closed  and  Variable-Volume  Simulators",  University  of 
Illinois  Teclinical  Report  AAE  73-6,  September  1973. 
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the  bomb  at  the  beginning  of  the  experiment.  Formulas  for  that  case 
are  given  in  Section  5.  These  formulas  show  how  the  pressure  rise  in 
the  bomb  is  influenced  by  the  burning  rate  and  by  the  experimental  set- 
up, respectively.  By  controlling  the  latter,  the  computation  of  the 
burning  rate  can  be  greatly  simplified. 

An  example  for  the  application  of  the  new  formulas  to  closed  bomb 
measurements  is  given  in  Section  6. 


2.  ASSUMPTIONS 

The  details  of  the  combustion  process  and  of  the  associated  gas 
flow  in  a closed  bomb  experiment  are  extremely  complicated.  For  the 
present  purpose  we  shall,  however,  neglect  particulars  of  that  flow  and 
assume  that  space  average  values  of  state  variables  describe  the  state 
of  the  gases  in  the  bomb  at  any  time  sufficiently  accurate.  We  shall 
also  neglect  the  kinetic  energy  \ersus  the  internal  energy  of  the  gas«s. 
These  assumptions  may  be  summarited  by  the  following  statement  (Assump- 
tion A) : 

(A)  Combustion  gases  mix  with  gases  present  in  the  bomb  in  a time 
which  is  short  compax'cd  with  the  combustion  time. 

Our  analysis  will  not  cover  any  chemical  reactions  and  wo  want  to 
describe  the  combustion  gas  by  a single  equation  of  state.  Wo  make, 
therefore,  the  following  two  assumptions: 

(B)  The  composition  of  the  combustion  gas  is  independent  of  the 
coiuiitions  in  the  bomb. 

(C)  Chemical  gas  phase  reactions  take  place  only  within  a region 
arbitrarily  close  to  the  surface  of  unbumt  propellant  and  the  gas 
released  from  that  region  has  a constant  (flame)  temperature. 

We  intend  to  consider  a single  homogeneous  propellant  only  for 
which  we  assume: 

(0)  The  density  and  the  temperature  of  the  unburnt  propel Jant  are 
constant. 

The  assumption  (D)  eliminates  from  our  considerations  possible  effects 
of  a preheating  of  the  propellant  during  the  test.  We  can  suaaarisc 
the  assumptions  (B),  (C)  and  (0)  as  follows: 

(h)  The  energy  released  by  the  combustion  is  proportional  to  the 
volume  consumed. 

We  assume  also: 

(F)  Heat  losses  to  the  walls  of  the  bomb  are  negligible. 
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This  assumption  is  probably  the  most  restrictive  one  in  our  analysis. 

Its  purpose  to  simplify  the  theory.  For  similar  reasons  we  require  a 
specific  form  of  the  equation  of  state,  namely: 

(G)  The  mixture  of  gases  in  the  bomb  obeys  an  equation  of  state 
of  the  form, 

p = E • f(m,V)  (1) 

where 

p is  pressure 

E is  the  internal  energy  of  the  gases 

m is  the  gaseous  mass 

V is  the  volume  available  for  the  gases,  i.e.,  the  inital 
void  volume  plus  the  volume  of  propellant  burned. 

The  next  assumption  is  a standard  assumption  about  the  burning  pro- 
cess . Let 

r be  the  "linear  burning  rate",  which  we  seek  to  determine, 

X be  the  linear  distance  burnt, 

S(x)  be  the  burning  surface  of  the  propellant,  and 

t bo  the  time. 

The  assumption  is  then  as  follows: 

(H)  The  combustion  of  the  propellant  obeys  the  equations 

dx  « rdt 
atul 

dV  ■ S(x)rdt. 

Equations (2)  express  the  increase  of  the  volume  V,  which  is  available 
for  gases  in  the  bomb  due  to  the  combustion  of  the  propellant,  the 
corresponding  change  of  the  utibumt  propellant's  volume  is  -dV. 

(I)  The  propellant's  surface  S(x)  is  constant. 

The  last  assumption  (I)  is  satisfied,  for  instance,  in  the  experi- 
ments reported  in  Reference  2.  It  is  essentially  a requirement  on  the 
geometry  of  the  propellant. 


S.  OUTLIHE  OF  TliE  KETHOD 

Based  on  the  assumptions  of  Section  2 we  may  describe  the 
combustion  process  in  the  closed  bomb  by  the  following  set  of  equations: 
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(3) 

(4) 

(5) 


dE  = CE  - p)  dV, 
dV  = Srdt, 

m = iEq  + (V  - Vo)Pp, 

A 

where  E is  a constant  representing  the  chemical  energy  released  per  unit 
volume  burnt,  is  the  mass  of  the  inert  gas  present  in  the  bomb  at 
ignition  time  (t  = o),  p is  the  density  of  the  unburnt  propellant,  and 
Vq  is  the  initial  volume^ available  for  the  gases.  Substitution  of  (5) 
into  (1)  yields  a relation  between  pressure,  energy  and  volume,  which 
is  of  the  form 


p = E • f (V) . 


(6) 


Substitution  of  (6)  into  (3)  yields 


dE  = (E  - E -f  (V))  dV 


(7) 


This  is  a linear  differential  equation  for  the  internal  energy  in  the 
closed  bomb.  Its  solution  is 

■V 


E(V)  = exp  [-F(V)] 


Eo  + B 


where  Eq  is  the  initial  value  of  E,  and 

.V 


i 


exp  [F(v)]dv 


F(VJ 


f (v)dv 


C8) 


(9) 


Substituting  (8)  into  (6J  we  obtain  a relation  between  the  pressure 
and  volume  in  a closed  bomb: 


p(V)  - f(V)  exp[-F(V)3  Po/f(Vo)  + E / exp[F(v)]dv  ,(10) 


where  p-  is  the  initial  pressure  in  the  bomb. 


This  relation  ;ls  valid  for  any  closed  bomb  combustion  process,  which 
satisfies  the  assumptions  (A)  through  (G) , Note,  that  the  assumptions  (H) 
and  (I)  about  the  burning  process  have  not  been  utilized  to  derive  Equation 
(10).  The  equation  of  state  enters  Equation  (10)  through  the  function  f(V), 

In  some  cases,  particularly  if  the  quadrature  (9)  cannot  be  carried 
out  analytically,  it  is  preferable  to  use  instead  of  Equation  (10)  the 
corresponding  differential  equation  for  the  pressure.  That  equation  can 
be  derived  as  follows.  Differentiating  (6)  we  obtain 


dp  = f(V)dE  ♦ Ef‘(V)dV. 


(11) 
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Next,  we  eliminate  E and  dE  between  (3),  (6),  and  (11),  rbtaming 


r 


j 

i 

i % 


I \ 


dp  = [f(V)(E-p)  + pf'(V)/f(V)]dV  (12) 

or 

p = [f(V)(g-p)  + pf(V)/f(V)]V  (13) 


(By  dots  we  indicate  the  derivatives  with  respect  lo  time.)  Equation 
(10)  is,  of  course,  the  integral  of  equation  (12).  Equation  (13)  may  be 
convenient,  if  the  burning  rate  formula  contains  the  time  derivative  p. 

Let  the  burning  rate  be  expressed  by  a function  containing  p,  p and 
a parameter  vector  y,  say, 

r = G(p,p,y)  (14) 


Then,  according  to  the  assumptions  (H)  and  (I) , we  may  express  the 
volume  by  the  formula  t 


V(t) 


(15) 


If  p(t)  has  been  observed,  we  can  rompute  numerically  V(t)  by  (IS) 
for  any  given  set  of  parameters  y.  Therefore,  the  right  hand  .sides  of 
equations  (10)  and  (13)  can  be  computed,  too,  as  functions  of  t for  any 
set  y.  Since  al;  the  left  hand  sides,  p(t)  or  ]ji(t),  of  those  equations 
have  been  observ:  u,  can  be  computed  from  observations,  we  obtain  a 

curve  fitting  problem.  The  constraint  is  thereby  either  Equation  (10) 
or  Bqaution  (13)  and  the  fitting  parameters  are  the  components  of  y. 


This  conceptually  simple  process  can  be  easily  extended  to  equations 
of  state  which  are  more  general  than  Equation  (1)  and  to  cases  where 
assumptions  (F)  and  (I)  do  not  hold.  However,  even  under  the  present 
restrictions  the  computations  might  require  non-trival  numerical  processing, 

such  as  extensive  tabulation  of  the  functions  F(V)  and  Jrdt  for  different 

parameter  values  and  repeated  numerical  substitutions.  The  algorithm  is 
simpler  if  some  of  the  quadratures  in  Equations  (9),  (10),  and  (1$)  can 
bo  carried  out  analytically.  The  corresponding  formulas  are  more  convenient 
for  discussion.  Wo  shall,  therefore,  derive  such  formulas  for  the  case 
of  NoblC“Abel  equation  of  state  with  constant  specific  heats.  Since 
that  equation  is  widely  used  in  interior  ballistics,  those  formulas  have 
also  their  own  merit  as  working  formulas. 


4.  FORMULAS  FOR  NOBLE-ABEL  GASES 

The  combustion  bomb  contains  according  to  our  assumptions  two  gases, 
namely,  an  inert  gas  at  the  beginning  of  the  experiment  to  which  the 
combustion  gas  is  added  by  tlie  burning  process.  Let  both  gases  obey 
Noble*Abei  equations  of  state  and  have  constant  specific  heats.  They 
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can  then  both  be  characterized  by  the  following  set  of  parameters; 

M - the  molar  mass, 

Y - the  ratio  of  specific  heats, 
n - the  co-volume, 
m - the  gas  mass. 


We  attach  the  index  ’’zero"  to  these  quantities,  if  they  refer  to 
initial  conditions,  that  is,  to  the  inert  gas,  and  the  index  "one”,  if 
they  refer  to  the  combustion  gases.  Thus,  m^  is  the  mass  of  gas  in  the 
bomb  at  ignition,  and  m^Ct)  (zero  at  ignition)  is  the  mass  of  combustion 
gases.  The  total  gaseous  mass  in  the  bomb  is 


m(t)  * m^  + m^Ct). 


(16) 


The  equations  of  state  for  the  mixture  of  both  gases  are  the  Noble-Abel 
equation 


ra  /Mj  + m /M 


(17) 


and  the  energy  equation,  which  is,  in  case  of  constant  specific  heats, 


E = RT 


. 1 , , ^ , - 1 \ 

«o  ifo  - V ■ 


(18) 


where  R is  the  universal  gas  constant  and  T is  the  temperature  of  the  mix- 
ture. Eliminating  RT  between  equations  (17)  and  (18),  we  obtain 


p = E 


o,/M,  + n 
II  o o 


vi-Vo) 


(19) 


We  note,  that  equation  (Id)  is  of  the  required  form  (1). 


At  this  point  it  is  convenient  to  introduce  instead  of  the  volume 
V a new  variable  z,  which  represents  the  volume  fraction  burnt.  Let  z 
be  defined  by 


z .8  (V  - V^)/Vp,  (20) 

where  is  the  initial  void  volume  and  V is  the  initial  volume  of  the 
unburnt  propellant.  As  the  combustion  progresses,  z varies  aonotonically 
from  zero  to  one.  In  terms  of  z we  can  express  the  void  volume  of  the 
chamber  by 


V “ tV„  ♦ 

P o 


(21) 
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i 

\ 


1 


t 


and  its  differential  by 

dV  » Vpdz . (22) 

Th^  gaseous  mass  m.  released  by  the  combustion  (or  the  solid  mass  consumed 
by  the  combustion)  is  in  terms  of  z given  by 

m^  « zffip  , (23) 

where  m is  the  initial  mass  of  propellant  placed  in  the  bomb. 

r 

Substituting  V and  m,  from  Equations  (21)  and  (23)  into  Equation  (19) 
we  obtain  after  simple  algebra 


“o  . 

* * ■_  ST 


L “o  . "i . yiW  ci-Vp)Vp~\ 

V ^ ^ 7 


(24) 


Let 


a » 


and 


(V^ 

m/itr-l) 

' " Wo^ 


- z ♦ b 

^ uVcr'(S7iT  • 


Vo 


(25) 


Equation  (24)  can  then  be  expressed  by 
p o E • h(z) 

P 


(26) 


(27) 


The  parameters  b anc  c of  the  function  h(z)  aro  positive  a:id  generally 
small  compared  to  one,  because  usuolly 


B <<  m 
0 “p 


(28) 


The  parameter  d can  be  negative  as  well  as  positive,  depending  on  the 
sign  of  1 - HjPp.  The  only  limitation  for  d is 

d > -1.  (29) 


4 
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This  limitation  is  a consequence  of  our  assumption  that  the  gas  mixture 
obeys  Noble-Abel's  equation  of  state  (19).  That  assumption  implies 
that  the  denominator  in  Eq.  (19)  remains  positive  during  the  complete 
combustion.  Therefore,  we  have  at  ignition  the  relation: 

V - n m„  > 0 (30) 

0 0 0 

and  after  completed  combustion  the  relation 

V + V - n m - n,m  > 0.  (31) 

o p 0 0 Ip 

Combination  of  these  inequalities  yields  the  inequality  (29) . 


Substituting  the  expressions  (22)  and  (27)  into  Eq,  (3)  we  obtain 

dE  = [EVp  - Eh(z)]dz.  (32) 

Eq.  (32)  corresponds  to  Eq.  (7),  where  the  absolute  volume  V instead  of 
the  relative  volume  z was  used.  With  h(z)  given  by  Eq.  (26)  the 
differential  equation  (32)  can  be  solved  analytically  except  for  one 
quadrature. 


The  initial  value  for  the  solution  of  the  differential  equation  (32) 
follows  from  Eq.  (27): 


E(0)  = Eq  = Pq  h(o) 


V c 
■ o ab 


^0  e 


1-Vo 


(33) 


For  the  analytic  solutions  of  Eq.  (32)  we  have  to  consider  the 
cases  d = 0 and  d ° 1/c  separately.  In  the  general  case  d 0 and 
d / 1/c  we  define 


I (34) 


A = u 


1-cd  ’ 


R = i 

® “ d TTd 


The  energy  E(z)  is  then  given  by 


E(z) 


S 


1 

(z+c)^(zd4l)® 


(35) 
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The  corresponding  pressure  p(z)  is  with  Eqs.  (26)  and  (27) 


C/--A  a z + b 
P(  ) “ ^ '^p  (z+c)  (zd+1) 


(z-i-b)a 

(z+c)^'^\zd+l)®*^ 


Po  * E 


(s+c)^(sd+l)®d 


Eq.  (36)  permits  to  compute  the  pressure  in  the  chamber  as  a function  of 
the  relative  volume  consumed.  It  depends  only  on  the  initial  pressure 
and  on  the  gas  parameters  listed  at  the  beginning  of  this  Section.  It 
is  independent  of  the  geometry  of  the  propellant  and  of  the  burning  rate 
equations . 

The  chemical  energy  constant  E is  equal  to  the  internal  energy  of 
the  combustion  gas  at  flame  tanperature,  per  unit  volume  of  the  unburnt 
propellant.  It  is  for  a Noble-Abel  gas  equal  to 


“'‘‘flane 


“p " Tpr  "p  • 


where  f is  the  "force  constant"  of  the  propellant,  (Ref, 5,  p.  96),  expres- 
sed in  metres,  and  g is  the  standard  acceleration.  The  dimension  of  £ 
is  that  of  a pressure. 

In  order  to  derive  an  equation  corresponding  to  Eq.  (13)  we  follow 
the  steps  outlined  in  Section  3,  obtaining  first  the  equation 

p = [(E-p)h(z)  ♦ ph‘(z)  /h(z)]z.  (38) 

With  h(z)  from  Eq.(26)  we  compute  the  derivative 


f 2+V 

h*(z)  = 


(z*c)^(zd^l)^ 


and  substitute  h(z)  and  h'(z)  into  Eq.(38).  The  result  is 

A „ ffg.p)  .p  (40) 

P (z+cHztiil)  P (z^b)(z^c)(zd»l) 

Eq.(36)  gives  the  pressure  p as  a function  of  the  relative  volume 
consum^.  In  cases  of  practical  interest  p(z)  is  strictly  monotonic 
function  and  can  be  inverted  to  yield  z as  a function  of  p.  By  substitu- 
tion of  that  inverted  function  into  Bq.(40)  one  obtains  an  expression  of 
the  type 


Corner,  F.,  "Theory  of  the  Interior  Ballistics  of  Guns"  John  Wiley  aitd 
Sons,  Now  York,  19S0. 
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where  C(p)  is  independent  of  the  burning  rate  and  of  the  burning  surface 
conditions.  This  equation  is  a central  equation  for  the  nmerical 
evaluation  of  closed  bomb  measurements.  Usually  it  is  assumed  that  C(p) 
can  be  approximated  by  a constant,  namely,  the  final  pressure  in  the  bomb, 
Pend  l"^)*  The  equation 

P ■ P,„<i  • i (42) 

can  be  derived  theoretically  if  one  assumes  uniform  gas,  constant 
temperature,  negligible  volume  change,  Po‘^*«  P _j>  and  n,P„  ® 1. 

Another  approximation,  ® P 

P=Pend^*^'‘ 

where  b'  is  an  experimental  fitting  parameter,  is  used  in  Refernces  2 and 

3. 


In  our  analysis  C(p)  is  controlled  by  the  setup  of  the  experiment. 
Because  a constant  CCp)  greatly  facilitates  the  processing  of  the  data, 
it  is  desirable  to  plan  the  experiments  such,  that  a constant  is  a good 
approximation  to  C(p).  In  the  next  Section  we  shall  show  by  some  examples 
how  C(p)  depends  on  various  parameters  of  the  experiment. 

The  fitting  of  experimental  data  with  the  aid  of  Eq,(41)  can  be 
best  illustrated  by  an  example.  Let  the  burning  rate  be  given  by  the 
equation 


r o 


(44) 


whore  0 and  n arc  free  parameters  and  a is  the  thermal  diffusivity  of  the 
unburnt  solid  propellant.  The  first  term,  Bp”,  on  the  right  hand  side 
of  Eq,  (44)  is  the  usual  form  of  a "steady  burning  rate  equation".  The 
benefits  of  adding  the  second  ("transient") term  are  discussed  at  length 
in  Reference  3.  For  some  propellants  such  a term  results  in  a more 
realistic  description  of  the  burning  behavior,  particularly  for  low 
pressures.  The  time  derivative  i is  given  with  Eqs.  (15)  and  (20)  in  this 
case  by 


Substitution  of  this  expression  into  Eq.  (41)  yields  a differential 
equation  for  p(t),  namely, 


p » C(p)  ;^[6p”  + I * P 
P L P J 


(46) 


Obviri,’y  Cq.  (46)  can  be  integrated  analytically  if  C(p)  is  constant. 
The  i 'gral  is  an  equation  containing  p,t  and  the  unknown  parameters 
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3 and  n.  That  equation  can  be  used  as  a constraint  functional  for  fitting 
the  observed  p(t). 

For  sake  of  completeness  we  note  here  also  the  final  formulas  for 
the  special  cases  excluded  above. 


First  we  consider  the  case  d=0,  that  is»n,p  = 1.  The  pressure 
function  p(z)  is  in  this  case  


PCz) 


(z*b)ae 


-az 


(z+c) 


(b-c)a+l 


, Cb“c)a+1 
ab 


. E J (s*c) 


The  relation  between  p and  z,  corresponding  to  Eq.  (40),  is 
p = [Ea  + (c-b-a)p3  ^ i. 

In  the  special  case  d=l/c  we  obtain  the  corresponding  equations 


(48) 


P(z) 


and 


(z-t-bjac^*^^  / (c-b)az 


(47) 


(^)[^  * ^ if 


(e"-p)  . p 


(Z+C) 


(z+bj  (z-^c)' 


z . 


(50) 


S.  DISCUSSION  OF  A SPECIAL  CASE 

The  formulas  derived  in  the  previous  section  still  contain  integrals 
which  must  be  evaluated  numerically  in  order  to  compute  the  function  C(p), 

In  this  section  we  shall  discuss  a special  case  which  produces  simplifi- 
cations of  those formulas.  The  purpose  of  the  discussion  is  to  provide 
a general  picture  of  the  properties  of  C(p).  Particularly  we  are  interested 
to  find  conditions  under  which  C(p)  can  be  approximate!  by  a constant. 

A common  assumption  for  the  treatment  of  closed  bomb  problems  is  to 
assume  ideal  gases  (no‘nt °0) . In  our  analysis  such  an  assumption  does 
not  simplify  any  expressions. 

Inspecting  the  equations  of  the  previous  section  we  observe  that 
considerably  simpler  formulas  can  be  obtained  if  the  parameters  b and 
c of  the  function  h(z)  can  either  be  neglected  or  are  equal.  The  latter 
is  the  case,  for  instance,  if  the  inert  gas  and  the  combustion  gas  are 
identical.  However,  in  order  to  make  b»c,  it  is  sufficient  to  require 
only 

, Yq  “ Y|-  (51) 
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We  assme  for  the  rest  of  this  Section,  that  Eq.(51)  holds.  The  function 
h(s)  is  then  simply  given  by 


h(z)  = 


(52) 


and  the  quadratures  of  Section  4 can  be  carried  out  analytically.  Thereby 
we  have  to  consider  the  special  cases  n^Pp  = 1 and  tij^pp  = Yj^  separately. 

In  the  general  case,  that  is,  for  nj^P„  / 1»  Vj^p  / Yp  and 
Yo  = Yi  we  obtain  the  following  formulas  — 2 


pw  = = slr*[i’o-^5lr]  = 

- p * r F -§_l 

“ ^ B+l  [Pend  ’ ^ B+1  J \zd+l / ’ 


(53) 


(B+2)/(B+l) 

l(B+2)/(B*l) 


(54) 


In  the  special  case  n^^Pp  ° Y^  and  y^  ° Y^  we  have 
P(5i)  = Pp  - E • In  (zd+1)  B 


In  the  case  n^^Pp  “ 1 and  Yq  “ Yj  the  formulas  are 
p B E - (E  - Pp)  exp  (-az)  » 

- i ’ ?end^  ®*P 


(55) 


(56) 


(57) 
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and 


C(p)  = a • (E  - p) . 


(58) 


Typical  examples  of  the  functions  p(z),  C(p)  and  of  the  temperature 
within  the  bomb  are  shown  in  Figures  1 through  3.  They  illustrate  the 
dependence  of  these  functions  on  some  parameters  of  the  closed  bomb 
experiment.  The  basic  parameter  values  used  for  these  calculations  are 
listed  in  Table  1. 


The  amount  of  propellant  enters  the  analysis  in  form  of  the  volume 
ratio  V /V^  (see  Eq.(25)).  A more  convenient  measure  for  the  amount  of 
propellant  is  the  ratio  of  the  propellant's  volume  to  the  volume  of  the 
bomb.  That  ratio  is  chosen  as  the  variable  parameter  for  the  curves  in 
Figure  1.  It  is  related  to  V /V^  by  the  equation 


1-V^/V, 


p bomb 


(59) 


The  examples  shown  in  Figure  1 indicate,  that  an  approximation  of  C(p) 
by  a constant  is  feasible  in  the  present  case  only  for  loading  volume 


ratios  less  than  0.1.  For  these  cases  and  about  one  atmosphere 

(lOSpa)  initial  pressure  C*  = is  a reasonable  approximation  to  C(p). 
The  corresponding  temperature  curves  T(p),  are  almost  constant  over  most 
of  the  pressure  range,  if  Vp/V|j  >0.1.  For  high  loading  volume  ratios 
the  temperature  in  the  bomb' reaches  a maximum  at  lower  pressures  and  then 
decreases,  as  the  combustion  progresses. 


The  examples  in  Figure  2 show  the  dependence  of  the  functions  C(p) 
and  T(p)  on  the  inital  pressure  p^.  It  is  interesting  to  note  that  for 
small  loading  volume  ratios  (Vp/V{jgjj,^j  = 0.05  in  Figure  2)  the  function 
C(p)  can  be  best  approximated  by  a constant  which  is  less  than  Pond- 

Figure  3 shows  the  significance  of  the  co-volume  term  for  the  analysis. 
The  special  cases  “ l/p«  and  n^  “ v/pp  exhibit  no  special  features  in 
these  examples.  ^ ^ 


An  inspection  of  the  formulas  of  this  section  shows  that  in ^general 
C(p)  can  be  approximated  by  a constant  if  p <<  0.  The  value  of  E depends 
on  the  physical  parameters  of  the  propellant  only.  The  maximum  value  of 
the  pressure  depends  on  the  other  hand  also  on  the  ratio  Vp/Vg  as  well 
as  on  pp.  Therefore,  the  condition  p <<  § limits  for  a given  propellant 
the  maximum  admissible  ratio  Vp/V^  and  the  maximum  initial  pressure  Pp. 

The  examples  presented  indicate  also  that  the  behavior  of  the  tempera- 
ture in  the  bomb  is  by  itself  not  relevant  to  the  question  whether  C(p) 
can  bo  approximated  by  a constant. 


{ 

I 

I 
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Table  1 


Basic  Parameter  Values  for  Figures  1 Through  3 


Inert  gas 


Combustion  gas 


Yo  = 1.21 

Hq  = 1.375  • 10'^  m^/kg 

Mq  = 28  • lo"^  kg/mol 

T^  = 300  K 
o 

Pq  = lO^Pa 


Yi  = 1.21 

= 1.09  • 10’^  m^/kg 
f • g = 1.186  • 10^  m^/s 

Aflame  = 

p = 1.62  • 10^  kg/m^ 

P 

(E  = 9.15  • 10^  Pa) 


Figu'*  3.  Closad  Bomb  FuncMons  for  Difforoni 
Co’volumos 

Fiiod  Poromortrt;  V^/V«{)mi*0.3 
P,»IO*Po 

VoriobU  Poromoror.' 1),  •O;O.3’10r®il//if^iX.  //»^jtO**|m*/ko) 


6.  EXAMPLE 


In  this  Section  we  shall  give  an  example  f the  application  of  the 
closed  bomb  data  reduction  algorithm.  The  data  are  taken  from  tm  recent 
experiments  at  the  Ballistics  Research  Laboratories  (Ref.  6).  The  important 
parameters  of  the  experiments  are  listed  in  Table  2.  Because  in  this 
case  Yq  ® Yi»  the  simple  formulas  of  Section  5 could  be  used  to  compute 
the  function  C(p),  The  result  shown  in  Figure  4 indicates  that  C(p)  is 
within  10%  of  the  end  pressure  throughout  the  combustion  process.  However, 
the  end  pressure  computed  by  Eq.  (53)  was  bigger  than  the  observed  end 
pressure.  It  was  conjectured  that  the  reason  for  this  discrepancy  could 
be  an  overestimation  of  the  propellants  ’•force"  constant.  Assuming 
therefore  a 34.8%  smaller  "force"  constant  both  pressures  were  matched. 

The  corresponding  new  C(p) -curve  was  found  to  be  within  plotting  accuracy 
of  Figure  4 equal  to  the  curve  computed  for  the  higher  "force"  constant, 
although  the  corresponding  temperatures  in  the  bomb  differed  by  up  to 
1200K.  It  was  then  concluded  that,  independently  of  the  true  value  of 
the  "force"  constant,  the  observed  p^nj  is  a reasonable  approximation  to 
C(p). 


Let  C*  be  a constant  approx j,..at ion  to  C(p).  The  integrated  constraint 
equation  (46)  is  then  for  n ^ 1 


FCt.PiS.n)  = jEj.  . i . (p‘-“ 


a , -2n  -2n,  ^ ^ . 

+ ^ (P  -Po  3 - t + to  = 0, 
28 


(60) 


and  for  n = 1 


F(t,p;B,l)  = J (P/Po^  * 

+ (p'^-Po*^)  - t to  = 0 , 

2a 


(61) 


In  Eqs.  (60  and  (61)  to  and  Pq  is  a pair  of  corresponding  values  of  ob 
served  time  and  pressure,  respectively. 


In  the  present  case  the  terms  containin'?  the  thermal  diffusivity 
u were  found  to  be  insignificant  and  were  therefoienot  used  in  the 
final  evpluations.  Also,  the  following  more  "eneral  steady  state 
burning  rate  formula  was  used  for  the  analysis: 

r = B(P  - p)".  (62) 


^Knapton,  J.D.,  Interior  ballistics  Laboratory,  USA  Ballistic  Research 
Laboratories,  private  communication.  The  experiments  have  IBL  numbers 
167-3  and  167-4,  dated  21  March  1974.  The  propellant  sample  was  JPN, 
coated  wit.h  a polystyrene  surface  inhibitc  ,. 
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Table  2 

Parameters  of  the  Experiment 


V /V.  . = 0.1278 

o'  bomb 


C*  =» 


^end 


= 2.04  • 10®  Pa 


V /S  = 6.4922  • 10'^  m 
P 


Inert  gas 


To  = 1-21 

n^  s 1.37S  ‘ 10’^  m^/kg 
Mq  = 28  • lO’^  kg/mol 
Tq  « 300  K 
p^  « 1.013  • 10®  Pa 


Combustion  gas 

Yj  » 1.21 

B 1.093  • 10"®  m®/kg 
Mj  » 25.66  * 10"®  kg/mol 

p a 1.62  • 10^  kg/m® 

^ 6 2 2 
Given  fg  ■ 1.186  * 10  m /s 

Used  fg  a 0.7738  • 10®  m^/s^ 

(E  « 5.970  • 10®  Pa) 


29 


With  this  burning  rate  formula  the  integrated  constraint  functional  is 
for  n / 1 

V 


FCt,p;3,n,p)  = ^ I Cp-p)*”'  - CPo“P)*'“  - t + t^  = 0 (63) 

and  for  n = 1 


,1-ii 


,1-n 


(64) 


The  equations  (63)  and  (64)  define  a three  parameter  family  of  curves 
in  the  t,p-plane.  They  all  pass  through  the  node  (tojpg),  which  corres- 
ponds to  an  arbitrary  integration  constant.  That  constant  is  an  additional 
parameter  of  the  curves  and  must  be  included  in  the  curve  fitting  process. 
In  the  present  case  it  is  convenient  to  chose  an  arbitrary  Pq  and  treat 
tg  as  a free  par^eter.  The  adjustment  problem  has  then  a total  of  four 
parameters:  6,n,p,  and  t^. 


The  pressure  data  were  obtained  simultaneously  from  two  pressure 
transducers  (Figure  5).  Consequently  ter  each  time  value  two  pressure 
values  were  recorded.  These  double  recordings  were  used  to  assess  the 
accuracy  of  the  measurements.  By  inspection  of  the  data  scattering 
and  of  the  discrepancies  between  simultaneous  pressure  readings  it 
was  estimated  that  the  standard  error  of  each  pressure  reading  was 
approximately 

6p  = 2 • 10^  + 0.02p  [Pa].  (65) 

The  standard  errors  of  the  time  observations  were  assumed  to  be 


e^  e 0.1  [ms].  (66) 

Correlations  between  pressure  and  time  readings  were  assumed  to  be 
negligible. 

For  the  data  reduction  example  pressure  readings  were  sampled  from 
a more  detailed  list  at  O.Sms  or  1.0ms  intervals.  Figure  6 shows  the 
sampled  pressure-timo  data  from  two  experiments  together  with  some 
3e-error  ellipses,  Indicating  the  accuracy  of  the  data.  High  pressure 
measurements  appear  in  Figure  6 to  be  more  accurate  because  the  pressure 
is  plotted  in  logarithmic  scale.  In  fact  only  the  relative  accuracy 
is  better  at  high  pressures. 

The  rapid  pressure  Increase  shown  in  Figure  6 occured  in  the  experi- 
ments about  1.5s  after  ignition.  Pressure  readings  prior  to  that  Increase 
have  a very  high  noise  level  and  cannot  be  used  without  special  treatment. 
In  the  present  evaluation  examples  we  did  not  use  those  observations  at 
all.  The  complete  p(t)  curve  has  an  appearance  which  is  schematically 
shown  in  Figure  7.  Only  a small  part  of  the  curve  for  t l>5s  is  well 
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defined  by  measurements.  The  time  coordinate  of  the  ignition  point  is 
known  only  approximately  and  data  between  t * 0 and  t = 1.5  are  practically 
not  available. 

The  constant  part  of  the  curve  (after  completed  combustion)  is  not 
described  by  eqs.  (63)  and  (64).  The  corresponding  data  should,  therefore, 
not  be  used  for  data  reduction  (except  for  determining  Pend) • ^n  the 
other  hand  the  transition  from  slow  to  fast  pressure  increase  at  about 
1.5  s is  not  modeled  either  by  the  family  of  curves  defined  by  eqs.  (63) 
and  (64) . It  seems  that  in  the  experiment  two  different  burning  processes 
take  place:  one  for  pressures,  say,  below  20MPa  and  a different  one  for 
higher  pressures.  Under  these  circumstances  low  pressure  observations 
and  high  pressure  observations  should  be  evaluated  separately.  We  have 
therefore  excluded  from  our  analysis  those  observations  which  are  obviously 
within  the  slow  pressure  increase  region  or  within  the  transition  between 
slow  and  fast  increase  regions.  The  choice  was  made  by  visual  inspection 
of  Figure  6. 

The  restriction  of  data  to  the  interval  between  about  20MPa  and  200MPa 
means  that  any  burning  laws  computed  from  those  data  are  valid  within  that 
pressure  region  only. 

Examples  of  curve  fitting  using  eqs,  (63)  and  (64)  are  shown  in  Figures 
8 and  9.  The  adjustment  was  obtained  by  the  non-linear  least  squares 
algorithm  described  in  Reference  7.  The  data  points  which  were  used  for 
the  curve  fitting  are  marked  in  Figures  8 and  9 by  black  squares.  From 
each  experiment  58  such  points  were  sampled  at  0.5  ms  time  intervals 
within  the  fast  pressure  increase  region.  The  figures  show  the  fitted 
curves  together  with  their  three  standard  error  confidence  limits  for 
different  burning  rate  formulas,  namely, 

r o Bp,  r = Bp”,  r = B(p  - p) , and  r = 8(p  - p)”. 

It  appears  from  the  plots  that  the  curve  fitting  does  not  improve  if  two  - 
and  three-parameter  burning  rate  functions  are  used  instead  of  the  simple 
linear  function  r = Bp.  The  numerical  results  are  summarized  in  Table  3 
and  confirm  that  impros.sion.  The  values  of  W (the  weighted  sum  of  correction 
squares)  and  (tlte  square  root  of  variance  of  weight  one)  are  affected 
only  insignificantly  by  the  addition  of  parameters. 

The  reasons  for  this  result  become  clear  if  one  considers  the  magnitude 
of  the  correlation  coefficients  between  the  parameters.  They  are  all  very 
close  to  one,  particularly  the  correlation  coefficients  between  0 and  n. 

This  indicates  that  ail  parameters  affect  the  fitting  curve  in  a similar 
manner  within  the  range  of  observations. 
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Tattle  S.  teasi  &paxcs  Values  and  Standard  Errors  of  Parameters  in  the  timing  Rate  Formula  r = 0(p  - p) 


In  Figure  10  we  show  the  least  squares  values  of  B and  n in  a n,lgB 
plane.  The  strong  correlation  between  these  parameters  let  the  error 
ellipses  look  like  inclined  error  bars.  The  JPN  parameter  values  of  Ref. 

8 are  plotted  in  Figure  10  for  comparison. 

The  differences  between  burning  rates  deputed  using  the  various 
burning  rate  formulas  are  not  as  large  as  it  appears  considering  the 
values  of  the  parameters  only.  If  the  burning  rate  is  plotted  as  a 
function  of  pressure,  then  all  burning  rate  formulas  produce  comparable 
burning  rate  values  withing  the  range  of  observations.  (See  Figures 
11a  through  lid) . The  main  differences  between  the  various  burning  laws 
are  their  different  standard  errors.  The  errors  are  larger,  for  formulas 
with  more  parameters,  because  the  simplest  curve  (r  = gp)  is  already 
located  well  within  the  data  scattering. 

The  similarity  of  the  burning  laws  is  demonstrated  also  in  Figure 
12,  where  the  various  burning  rate  functions  of  Table  3 are  plotted  in 
one  graph.  The  eight  curves  corresponding  to  the  first  eight  burning 
rate  functions  of  Table  3 are  close  to  each  other  within  the  pressure 
interval  of  observations.  The  burning  rate  function  of  Reference  8 differs 
from  the  other  functions  by  a factor  of  25  to  50  within  the  same  pressure 
interval.  The  reason  for  this  difference  is  possibly  the  fact  tirnt 
Reference  8 deals  with  burning  rates  at  low  pressures. 

In  summary*  these  results  ihdicate  that  the  simplest  burning  rate 
formula  r =0p  is  a good  approximation  for  JPN  at  pressures  between 
20  MPa  and  200MPa.  The  coefficient  6 can  be  computed  as  a weighted 
average  of  the  corresponding  values  of  both  experiments.  The  final 
average  burning  rate  formula  is  then 

r o (3.431  ♦ 0.029)  10‘®  p,  (65) 

where  r is  expressed  in  m/s  and  p in  Pa. 


7.  CONCLUSIONS  AND  RBCOWENOATIONS 

The  closed  bomb  theory  presented  in  this  report  permits  one  to 
establish  rational  limits  for  the  validity  of  simplified  data  reduction 
formulas  which  are  presently  in  use.  Expansion  of  the  theory  to  more 
general  cases,  including  for  instance  variable  surface  propellants  and 
variable  C(p)  is  conceptually  simple,  and  requires  merely  the  development 
of  complicated  computer  codes. 

The  data  reduction  process  has  been  fomilated  as  a nonlinear 
curve  fitting  task.  Such  an  approach  permits  one  to  apply  standard  curve- 
fitting theory  and  algorittns  to  the  problem, whereby  valuable  information 
about  the  accuracies  and  correlations  of  the  burning  rate  parameters  can 
be  obtained.  It  was  found  that  in  general  these  parameters  are  strongly 
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correlated.  The  .standard  error  of  the  burning  rate  formula  was  found 
to  increase  rapidly  outside  the  observed  pressure  range.  Closed-bomb 
measurements  should  be  planned  therefore  so  that  a large  pressure 
range  is  covered  by  the  experiment,  including  particularly  those  pressures 
for  which  the  burning  rate  formula  is  to  be  used. 

The  burning  rate  of  JPN  was  determined  from  the  experiments  167-3 
and  167-4.  It  was  found  that  for  pressures  between  20MPa  and  200MPa, 

JPN  burning  rate  can  be  approximated  by  the  linear  function  r = 0p  with 
an  accuracy  compatible  with  the  accuracy  of  the  measurements.  The  burning 
rate  outside  the  above  mentioned  pressure  interval  cannot  be  determined 
from  the  given  data. 

The  least  squares  values  of  the  two  parameters  in  the  burning  rate 
formula  r = 3p^  have  a correlation  coefficient  very  close  to  minus  one. 
This  strong  correlation  indicates  that  one  of  the  ^rameters  is  redundant 
for  the  data  fitting  process  and  might  be  chosen  within  limits  by  other 
considerations.  In  the  present  case  n = 1 is  suggested  for  simplicity. 
Equally  good  agreements  with  data  could  be  obtained,  e.g.,  by  chosing 
n = 0.9orn  = 1.1. 

Several  general  improvements  of  the  described  closed  bomb  data 
reduction  process  are  possible.  Thus,  instead  of  evaluating  each 
experiment  separately,  all  experiments  involving  the  same  propellant 
could  be  evaluated  simultaneously.  The  least  squares  utility  programs 
presently  available  at  BRL  can  be  used  for  such  purpose.  The  processing 
requires  only  the  coding  of  some  new  subroutines  representing  constraint 
functionals . 

Methods  for  the  evaluation  of  the  noisy  measurements  and  for  a joint 
evaluation  of  data  for  all  pressures  can  also  be  developed.  In  connection 
with  such  a joint  evaluation  more  flexible  two-  and  three-parameter 
burning  rate  formulas  should  be  tested  numerically. 
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